Maximal width of the separatrix chaotic layer 



S.M. Soskin 1 ' 2 ' 3 and R. Mannella 4 

Institute of Semiconductor Physics, National Academy of Sciences of Ukraine, 03028 Kiev, Ukraine 

2 Abdus Salam ICTP, 34100 Trieste, Italy 
3 Physics Department, Lancaster University, Lancaster LAI 4YB, UK 
* Dipartimento di Fisica, Universita di Pisa, 56127 Pisa, Italy 

The main goal of the paper is to find the absolute maximum of the width of the separatrix 
chaotic layer as function of the frequency of the time-periodic perturbation of a one-dimensional 
Hamiltonian system possessing a separatrix, which is one of the major unsolved problems in the 
theory of separatrix chaos. For a given small amplitude of the perturbation, the width is shown to 
possess sharp peaks in the range from logarithmically small to moderate frequencies. These peaks 
are universal, being the consequence of the involvement of the nonlinear resonance dynamics into 
the separatrix chaotic motion. Developing further the approach introduced in the recent paper by 
Soskin ct al. [PRE 77, 036221 (2008)), we derive leading-order asymptotic expressions for the shape 
of the low-frequency peaks. The maxima of the peaks, including in particular the absolute maximum 
of the width, are proportional to the perturbation amplitude times either a logarithmically large 
factor or a numerical, still typically large, factor, depending on the type of system. Thus, our theory 
predicts that the maximal width of the chaotic layer may be much larger than that predicted by 
former theories. The theory is verified in simulations. An application to the facilitation of global 
chaos onset is discussed. 
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I. INTRODUCTION 

Separatrix chaotic layers (SCLs) play a fundamental 
role for Hamiltonian chaos and may be important in a 
broad variety of subjects in physics and astronomy [J, 0, 
i, 0, i, 0, 0, 0, 0, M, E3- 0ne of the most important 
characteristics of the layer is its width in energy or in 
related quantities. It can be easily found numerically by 
means of integration of the Hamiltonian equations with 
a set of initial conditions in the vicinity of the separatrix. 
But it is important also to be able to find it theoretically. 
There is a long and rich history of the corresponding 
studies. The results may be classified as follows. 

1. Heuristic analytic results. 
Consider a ID Hamiltonian system perturbed by a weak 
time-periodic perturbation: 

H = H (p,q) + hV(p,q,t), (1) 
V(p,q,t + 2ir/u f ) = V(p t q,t), h < 1, 

where Hq (p, q) possesses a separatrix and, for the sake of 
notation compactness, all relevant parameters of Hq and 
V, except possibly Uf, are assumed to be ~ 1. 

There were a few heuristic criteria set by physicists (see 
e.g. 0, 0, 0, 0, 0, 0|) which gave qualitatively similar 
results for the SCL width AE in terms of energy E = 
H Q (p, q): 

AE = AE(u f ) ~ uj f S, (2) 
S=h\e\, 

|e| ~ 1 for uif ~ 1, 

|e| oc exp(— au>f) <C 1 (a ~ 1) for u>f 1. 



The quantity S = h\e\ is called the separatrix split Q 
(see also Eq. (4) below): it determines the maximum 
distance between the perturbed incoming and outgoing 
separatrices 0, 0, 0, 0, 0, 0, 0, 0, 0] . 

It follows from (2) that the maximum of AE lies in the 
frequency range ujj ~ 1 while the maximum itself is ~ h: 

AE max = max{AE(Lj f )}^h, oA™" ~ 1. (3) 

2. Mathematical and accurate physical results. 

Many papers studied the SCL by mathematical and ac- 
curate physical methods. 

For the range ojf 3> 1, there were many works study- 
ing the separatrix splitting (see the review Q and refer- 
ences therein) and the SCL width in terms of the normal 
coordinates (see the review |9| and references therein). 
Though quantities studied in these works typically differ 
from those studied by physicists [l], 0, 0, 0, 0, 0| , they 
implicitly confirm the main qualitative conclusion from 
the heuristic formula (2) in the high frequency range: if 
LOf 3> 1 the SCL width is exponentially small. 

There were also several works studying the SCL in 
the opposite (i.e . adiabatic) limit ujf — ► 0: see e.g. 
[i~2l EH EH El) El| and references therein. In the 
context of the SCL width, it is most important that 
AE(ujf — > 0) ~ h for most of the systems [I2,E1)E3|- For 
a particular class of systems, namely for ac-driven spa- 
tially periodic systems (e.g. the ac-driven pendulum), 
the width of the SCL part above the separatrix diverges 
in the adiabatic limit [H, [la ]: the divergence develops 
for Lj f < l/ln(l/fr). 

Finally, there is a qualitative estimation of the SCL 
width for the range uif ~ 1 within the Kolmogorov- 
Arnold-Moser (KAM) theory [9j while the quantitative 
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estimation within the KAM theory appears to be very 
difficult for this frequency range [17]] . ft follows from the 
results in d that the width in this range is of the order 
of the separatratrix split while the latter is of the order 
of h. 

Thus, from the above results, it could seem to follow 
that, for all systems except the ac-driven spatially pe- 
riodic systems, the maximum of the SCL width is ~ h 
and occurs in the range Wf ~ 1, quite in agreement with 
the heuristic result (3). Even for the ac-driven spatially 
periodic systems, this conclusion could seem to apply to 
the width of the SCL part below the separatrix, for the 
whole frequency range, and to the width of the SCL part 
above the separatrix, for uif ~ l/ln(l//i). 

3. Numerical evidences of high peaks in AE(ojf) 
and their rough estimates. 

The above conclusion does not agree with several nu- 
merical studies carried out du ring the last decade (see 
e.g. [H M, GJ, H M, SH SI SI) which have revealed 
the existence of sharp peaks in AE(ujf) in the frequency 

range l/ln(l//i) ~ uif ~ 1 the heights of which greatly 
exceed h (see also Figs. 2, 3, 5, 6 below). Thus, the 
peaks represent the general dominant feature of the func- 
tion AE(uif). The peaks were related by the authors of 
[ill QjJ S3, SIL S3, SI to the absorption of nonlinear reso- 
nances by the SCL. For some partial case, rough analytic 
estimates for the position and magnitude of the peaks 
were made in [l^, SI • 

4. Approach to an accurate description of the 
peaks. 

Accurate analytic estimates for the peaks were lacking. 
It is explicitly stated in the review [21j that the search 
for the mechanism of the involvement of resonances into 
the separatrix chaos and an accurate analytic descrip- 
tion of the peaks are being among the most important 
and challenging tasks in the study of separatrix chaos. 
The first step towards this accomplishment was done in 
the recent papers 0, SI , where a new approach to the 
theoretical treatment of the separatrix chaos for the rel- 
evant frequency range was developed and applied to the 
problem of the onset of global chaos between two close 
separatrices. An application of the approach to the, more 
common, single-separatrix cases was only discussed in 

[US!. 

The present paper formulates the basic ideas of the ap- 
proach in terms more general than [24], SI and, on the 
basis of this approach, develops the first ever accurate 
theoretical description of the peaks i.e. of the SCL width 
as a function of frequency in the range of the maximum 
of the width, which is the most important range from 
the physical point of view. In particular, we show that 
the maximal width of the separatrix chaotic layer may 
be much larger than it was assumed before. In the latter 
context, all systems are classified by us into two differ- 
ent types: for systems of type I, the ratio between the 
maximal width and the perturbation amplitude h loga- 
rithmically diverges in the asymptotic limit h — > while, 



for systems of type II, it asymptotically approaches a 
constant (still large, typically). 

Though the form of our treatment differs from typi- 
cal forms of mathematical theorems in this subject (cf. 
[1) II), the results yield the exact leading-order term in 
the asymptotic expansion of the width in the parameter 
of smallness a = l/ln(l//i). Our theory is in excellent 
agreement with the results of numerical integration of the 
equations of motion. 

Sec. II describes the basic ideas of the approach. Sec. 
Ill presents the classification into two types of systems, 
using rough estimates. Sec. IV develops the leading- 
order asymptotic theory for an archetypal example of 
type I and compares it with the numerical integration 
of Hamiltonian equations of motion. Sec. V develops 
the leading-order asymptotic theory for two archetypal 
examples of type II and compares it with the numerical 
integration. Next-order corrections are estimated in Sec. 
VI. Discussion of a few other issues, including in particu- 
lar an application to the global chaos onset, is presented 
in Sec. VII. Conclusions are drawn in Sec. VIII. 



II. BASIC IDEAS OF THE APPROACH 

The new approach, which is developed in [24L SI and 
here, may be briefly formulated as the matching between 
the discrete chaotic dynamics of the separatrix map in 
the immediate vicinity of the separatrix and the contin- 
uous regular-like dynamics of the resonance Hamiltonian 
beyond the close vicinity of the separatrix. The present 
section describes the general features of the approach in 
more details. 

The motion near the separatrix may be approximated 
by the separatrix map (SM) 0, SI, i, 0, 0, 0, 0, 0, [l8, 
SI, S3, SH, SI ■ It was introduced for the first time in [1 1 
and its various modifications were used in many studies 
afterwards, sometimes being called as the whisker map. 
It was re-derived in [26( rigorously, as the leading-order 
approximation of the motion near the separatrix in the 
asymptotic limit h — > 0, and an estimate of the errors 
was carried out too (see also [9| and references therein). 

We remind the main ideas which allow one to introduce 
the SM 0, SI, 0, H 0, i, 0, 1 Si SI, SI • For the sake of 
simplicity, let us consider a perturbation V that does not 
depend on the momentum: V = V(q,t). A system with 
an energy close to the separatrix value spends most of the 
time in the vicinity of the saddle(s), where the velocity 
is exponentially small. Differentiating E = Hq(jj, q) with 
respect to time and allowing for the equations of motion 
of the system (1), we can show that E = qdV/dq cx q. 
Thus, the perturbation can significantly change the en- 
ergy only when the velocity is not small i.e. during the 
relatively short intervals while the system is away from 
the saddle(s): these intervals correspond to pulses of ve- 
locity as a function of time. Consequently, it is possible 
to approximate the continuous Hamiltonian dynamics by 
a discrete dynamics which maps the energy E, the per- 
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turbation angle ip = uitt and the velocity sign a = sgn(q) 
from pulse to pulse. 

The actual form of the SM may vary, depending on 
the system under study, but its features, relevant in the 
present context, are similar for all systems. For the sake 
of clarity, let us consider the explicit case when the sepa- 
ratrix of Hq(p, q) possesses a single saddle and two sym- 
metric loops while V — qcos(ujft). Then the SM reads 



E. l+ i = E l + (Tihesm{tpi), 

Uf7r(3-sgn(E i+1 -E s )) 
<fi+i =<Pi-\ 7T~ TB — \ > 

a i+ i = Oi sgn(£; s - E i+ i), \a t \ = 1, 



(4) 



e = = 




SgI1 [ dp 




it—* — oc 


Ei= Ha(p,q)\ u _ 


Ifii = UJfti, 




Oi = sgn 




dp t 



dt 



dp 



sm(u> ft), 



A ' 



where E s is the separatrix energy, uj{E) is the frequency 
of oscillation with energy E in the unperturbed case (i.e. 
for h = 0), t{ is the instant corresponding to the i-th 
turning point in the trajectory q(t), and A is an arbi- 
trary value from the range of time intervals which greatly 
exceed the characteristic duration of the velocity pulse 
while being much smaller than the interval between the 
subsequent pulses [3, |, 1, 1, i, 1 0, ©, HI . 

Consider the two most general ideas of our approach. 

1. If a trajectory of the SM includes a state 
with E = E s and an arbitrary cp and a, then this 
trajectory is chaotic. Indeed, the angle ip of such a 
state is not correlated with the angle of the state at 
the previous step of the map, due to the divergence of 
uj^ 1 (E — > E s ). Therefore, the angle at the previous step 
may deviate from a multiple of 2tt by an arbitrary value 
and, hence, the energy of the state at the previous step 
may deviate from E s by an arbitrary value within the 
interval [— h\e\, h\e\]. The velocity sign a is not corre- 
lated with that at the previous step either (27] . Given 
that a regular trajectory of the SM cannot include a step 
where all three variables of the SM change random-like, 
we conclude that such a trajectory is chaotic. 

Though the above arguments appear to be obvious, 
they may not be considered as a mathematically rigorous 
proof, so that the statement about the chaotic nature of 
the SM trajectory which includes any state with E = E s 
should be considered as a conjecture supported by the 
above arguments and by the results of the numerical it- 
eration of the SM. Possibly, the mathematically rigorous 
proof should involve an analysis of the Lyapunov expo- 
nents for the SM (cf. [3|) but this appears to be a techni- 
cally difficult problem. We emphasize however that the 



rigorous proof of the conjecture is not crucial for the va- 
lidity of the main results of the present paper, namely of 
the leading terms in the asymptotic expressions describ- 
ing the peaks of the SCL width as a function of the per- 
turbation frequency. It will be obvious from the next item 
that in order to derive the leading term it is sufficient to 
know that the chaotic trajectory does visit areas of the 
phase space where the energy deviates from the separa- 
trix by values of the order of the ^pjarjitnxjsplit j5 = h\e\ 
which is a widely accepted fact 

2. As well known P,lli,i,[E[ 
at the leading-order approximation the frequency of 
eigenoscillation as function of the energy near the sepa- 
ratrix is proportional to the reciprocal of the logarithmic 
factor 




LO(E) 



\E-E.\ < AiJ = £ s 



b= 3-s g n(E-E^ (5) 



E s 



where E st is the energy of the stable states. 

Given that the argument of the logarithm is large in the 
relevant range of E, the function tv(E) is nearly constant 
for a substantial variation of the argument. Therefore, 
as the SM maps the state (Eq — E s ,(po,ao) onto 
the state with E = E\ = E s + a^he sin(<^o)? the value 
of oj(E) for the given sgn(croeshi(<^o)) is nearly the 
same for most of the angles ipo (except in the close 
vicinity of multiples of it), namely 



w(E)nu;?>, (6) 
ee u}(E s ± h), sgn(cr esin(v5 )) = ±1. 



Moreover, if the deviation of the SM trajectory from 
the separatrix increases further, u>{E) remains close to 
Wr^ provided the deviation is not too large, namely if 
ln{\E-E s \/h) < ln(AH/h). If uj f ~ U ±] , then the evo- 
lution of the SM (4) may be regular-like for a long time 
until the energy returns to the close vicinity of the sepa- 
ratrix, where the trajectory is chaotized. Such a behavior 
is especially pronounced if the perturbation frequency is 
close to u>r or uii ' or to one of their multiples of rela- 
tively low order: the resonance between the perturbation 
and the eigenoscillation gives rise to an accumulation of 
energy changes for many steps of the SM, which results 
in a deviation of E from E s that greatly exceeds the sep- 
aratrix split h\e\. Consider a state at the boundary of 
the SCL. The deviation of energy of such a state from 
E s depends on its position at the boundary. In turn, the 
maximum deviation is a function of u>f . The latter func- 
tion possesses the absolute maximum at ujf close to Lui +S> 

or u>r ' typically [28| , for the upper or lower boundary of 
the SCL respectively. This corresponds to the absorption 
of the, respectively upper and lower, lst-order nonlinear 
resonance by the SCL. 
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The above intuitive idea has been explicitly confirmed 
in [24]: it has been shown in the Appendix of [24j that, in 
the relevant range of energies, the separatrix map can be 
reduced to the system of two differential equations which 
are indentical to the equations of motion of the auxiliary 
resonance Hamiltonian which describes the resonance dy- 
namics in terms of the conventional canonically conjugate 
slow variables, action / and slow angle ip = nip — LOft 
where ip is the angle variable [3, H, [j, [H, H, 0] (see Eq. 
(16) below) while n is the relevant resonance number i.e. 

the integer number closest to the ratio Uf/ujj. . 

Thus, the result of the matching between the discrete 
chaotic dynamics of the SM and the continuous regular- 
like dynamics of the resonance Hamiltonian is the follow- 
ing [24]]. After the chaotic trajectory of the SM visits any 
state on the separatrix, the system transits in one step 
of the SM to a given upper or lower curve in the I — ip 
plane which has been labelled [24| respectivel y up per or 
lower generalized separatrix split (GSS) curve [29 ]: 

E = EW s $) = E a ±5\sm(j>)\, 6 = h\e\, (7) 

where S is the conventional separatrix split || while e is 
the amplitude of the Melnikov-like integral defined in Eq. 
(4) above (cf. 0, |, 1 1 i i 1 !, M, M M H S3), 
and the angle ip may take any value from one of the two 
ranges: either [Q,ir] or [it, 2tt] [3f3 |. 

After that, because of the closeness of uif to the n-th 
harmonic of ui(E) in the relevant range of E [3l|, for 
a relatively long time the system follows the nonlinear 
resonance (NR) dynamics (see Eq. (16) below), during 
which the deviation of the energy from the separatrix 
value grows, greatly exceeding 5 for most of the trajec- 
tory. As time goes on, ip is moving and, at some point, 
the deviation in energy from the separatrix value begins 
to decrease. This decrease lasts until the system hits the 
GSS curve, after which it returns to the separatrix just 
for one step of the separatrix map. At the separatrix, the 
slow angle ip is chaotized, so that a new stage of evolution 
similar to the one just described occurs, i.e. the nonlinear 
resonance dynamics starting from the GSS curve with a 
new (random) value of ip. 

Of course, the SM cannot describe the variation of the 
energy during the velocity pulses (i.e. in between instants 
relevant to the SM): in some cases this variation can be 
comparable with the change within the SM dynamics. 
This additional variation will be taken into account be- 
low, where relevant (see Sec. V below). 

One might argue that, even for the instants relevant 
to the SM, the SM describes the original Hamiltonian 
dynamics only approximately (26| and may therefore miss 
some fine details of the motion: for example, the above 
picture does not include small windows of stability on the 
very separatrix. However these fine details are irrelevant 
in the present context, in particular the relative portion 
of the windows of stability on the separatrix apparently 
vanishes in the asymptotic limit h — > 0. 



The boundary of the SM chaotic layer is formed by 
those parts of the SM chaotic trajectory which deviate 
from the separatrix more than others. As follows from 
the structure of the chaotic trajectory described above, 
the upper/lower boundary of the SM chaotic layer is 
formed in one of the two following ways [24], [25]]: 1) if 
there exists a self-intersecting resonance trajectory (in 
other words, the resonance separatrix) the lower/upper 
part of which (i.e. the part situated below/above the 
self-intersection) touches or intersects the upper/lower 
GSS curve while the upper/lower part does not, then the 
upper/lower boundary of the layer is formed by the up- 
per/lower part of this self-intersecting trajectory (Figs. 
1(a) and 1(b)); 2) otherwise the boundary is formed by 
the resonance trajectory tangent to the GSS curve (Fig. 
1(c)). It is shown below that, in both cases, the varia- 
tion of the energy along the resonance trajectory is larger 
than the separatrix split d by a logarithmically large fac- 
tor cx \n(l/h). Therefore, over the boundary of the SM 
chaotic layer the largest deviation of the energy from the 
separatrix value, AEsn) , may be taken, in the leading- 
order approximation, to be equal to the largest variation 
of the energy along the resonance trajectory forming the 
boundary, while the latter trajectory can be entirely de- 
scribed within the resonance Hamiltonian formalism. 

Finally, we mention in this section that, as obvious 
from the above description of the boundary, AEgtp} = 
Ai?!m (w/) possesses a local maximum ABitI : sm at ujf 
for which the resonance separatrix just touches the cor- 
responding GSS curve (see Fig. 1(a)). 

III. ROUGH ESTIMATES. CLASSIFICATION 
OF SYSTEMS. 

As obvious from Sec. II above, AE^^sm is equal, in 
the leading order, to the width AEnr of the nonlinear 
resonance which touches the separatrix. Let us make a 
rough estimate of AEnr'- it will turn out that it is possi- 
ble to classify all systems into two different types. With 
this aim, we expand the perturbation V into a Fourier 
series in t and a Fourier series in ip: 

V = ^^2V W (E,ip)exp(-Uv f t) + c.c. 
i 

= l^vW(E)exp(i(kip-lu; f t)) + c.c. (8) 

l,k 

As in the standard theory of a nonlinear resonance 

0, H, Hi S @> let us single out the relevant for a 

given peak, and denote its absolute value by Vb: 

V Q {E) = \V^\E)\. (9) 

Let us now roughly estimate the width of the reso- 
nance, using the pendulum approximation of the reso- 
nance dynamics @, i, i, H, H, 0] : 
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slow angle \j> 




slow angle \j> 




slow angle \j> 



FIG. 1: A schematic figure illustrating the formation of the 
peak of the function ABiJ(wj): (a) w/ = oj max ; (b) loj < 
^max; (c) ujf > o; max . The relevant (lower) GSS curve is 
shown by the dotted line. The relevant trajectories of the 
resonance Hamiltonian are shown by solid lines. The lower 
boundary of the layer is marked by a thick solid line: in (a) 
and (b) the lower boundary is formed by the lower part of the 
resonance separatrix while, in (c) it is formed by the resonance 
trajectory tangent to the GSS curve. The dashed line marks, 
for a given w/, the maximal deviation of the lower boundary 
from the separatrix energy E s . 



AE 



NR 



8hV Lu f 
\doj/dE\ ' 



(10) 



Of course, this approximation assumes the constancy 
of duj/dE in the resonance range of energies, while it 



is not so in our case: cu(E) oc l/\n(l/\E— EJ) in the 
vicinity of the separatrix [E H i, B, H 0, HI El H3, 
HI 0, so that the relevant derivative \du/dE\ ~ 
(uji^) 2 / (lu \E—E s \) strongly varies within the resonance 
range. However, for our rough estimate we may use the 
maximal value of \E — E s \, which is equal to AEnr ap- 
proximately. If ujf is of the order of w^' ~ ojo/ ^(l/h), 
then Eq. (10) reduces to the following rough asymptotic 
equation for AE^r'- 



AEnr 



V {E = E s ±AE NR )h\n(l/h), (11) 



The asymptotic solution of Eq. (11) essentially de- 
pends on Vq(E s ± AEnr) as a function of AEmr- In 
this context, all systems can be divided in the following 
two types. 

Type I. The separatrix of the unperturbed system has 
two or more saddles while the relevant Fourier coefficient 
y( L ) = y( L )(iJ,-0) possesses different values on adjacent 
saddles. Given that, for E — > E s , the system stays most 
of time near one of the saddles, the coefficient (E — > 
E s ,ip) as a function of ip is nearly a "square wave": it 
oscillates between the values at the different saddles. The 
relevant K is typically odd and, therefore, V~o(E — > E s ) 
approaches a well defined non-zero value. Substituting it 
in Eq. (11), we conclude that 



AE NR cx hln(l/h), 



0. 



(12) 



Type II. Either (i) the separatrix of the unperturbed 
system has a single saddle, or (ii) it has more than one 
saddle but the perturbation coefficient 7^' is identical 
for all saddles. Then V^ L \E — ► E s ,tfj), as a periodic 
function of significantly differs from its value at the 
saddle (s) only during a small part of the period in ip: this 
part is ~ uj{E)/u q - l/\n{l/\E s - E\). Hence, V (E S ± 
AE^r) oc \/ \q.{1/ AEnr). Substituting this value in Eq. 
(11), we conclude that 



AEnr oc h, 



0. 



(13) 



Thus, for systems of type I, the maximal width of the 
SM chaotic layer is proportional to h times a logarithmi- 
cally large factor oc \n(l/h) while, for systems of type II, 
it is proportional to h times a numerical factor. 

As shown below, the variation of energy in between 
the instants relevant to the SM is ~ h, which thus it is 
much less than AEmr (12) for the systems of type I and 
it is of the same order of AE^r (13) for the systems 
of type II. Therefore, one may expect that the maximal 
width of the layer for the original Hamiltonian system 
(1), AE^, is at least roughly approximated by that 

for the SM, AE^m , so that the above classification of 
systems is relevant to AE^' too. This is confirmed both 
by the numerical integration of equations of motion of the 
original Hamiltonian system and by the more accurate 
theory presented in the next two sections. 
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FIG. 2: Computer simulations for the ac driven pendulum 
(14) (an archetypal example of type I): the deviation AE^~^ 
of the lower boundary of the chaotic layer from the separatrix, 
normalized by the perturbation amplitude h, as a function 
of the perturbation frequency w/, for various h. The inset 
presents the same data but in logarithmic vertical scale and 
with the estimates by the heuristic Q, adiabatic [l3l ] and 
mathematical moderate-frequency 9] theories: the heuristic 
estimate is shown by the dotted line [32] while the adiabatic 
and moderate-frequency estimates are shown by the dashed 
line [33l ] . The inset explicitly shows that the simulation results 
exceed the estimates by the former theories by 1 or 2 orders 
of magnitude, for a wide range of frequencies. 



IV. ASYMPTOTIC THEORY FOR SYSTEMS 
OF TYPE I. 

For the sake of clarity, we consider a concrete example 
of type I, while the generalization is straightforward. 

Let us consider an archetypal example: the ac-driven 
pendulum (sometimes called as a pendulum subject to a 
dipole time-periodic perturbation) [J, [TH, : 



H = H Q + hV, 
P 2 

H = — - cos(g), 



V = —qcos(LUft), 



(14) 
/i < 1. 



Fig. 2 presents the results of computer simulations (i.e. 
of a numerical integration of the equations of motion) for 
a few values of h and several values olu)t. It shows that: 
1) the function AE(>(uif) indeed possesses sharp peaks; 
their height greatly exceed the estimates by the heuristic 
4 1 , adiabatic [1 31 ] and mathematical moderate- frequency 
9[ theories (see the inset); 2) as predicted by the rough 
estimates in Sec. Ill, the 1st peak of AE^ (uj f) shifts to 
smaller values of w/ while its magnitude grows, as h de- 
creases. Below, we develop the leading-order asymptotic 
theory and compare it with results of the simulations. 

Before moving on, we note that the SM (approximated 
in the relevant case by the nonlinear resonance dynam- 
ics) considers states of the system only at discrete in- 




0.2 O.i 0.4 0.5 0.6 0.7 0.8 0.9 





FIG. 3: An archetypal example of type I: ac-driven pendu- 
lum (14). Comparison of theory (solid lines) and simulations 
(circles): (a) the deviation AE^~\u) /) of the lower boundary 
of the chaotic layer from the separatrix, normalized by the 
perturbation amplitude h, as a function of the perturbation 
frequency ujj, for h — 10~ 6 ; the theory is by Eqs. (26), (31), 
(32), (38), (39) and (41). (b) the frequency of the 1st maxi- 
mum m as a function of h\ the theory is by Eq. 
(26). (c) the 1st maximum in AE^-~' (ojf)/h as a function of 
h; the theory is by Eqs. (34) and (26). 



stants. Apart from the variation of energy within the 
SM dynamics, the variation of energy in the Hamilto- 
nian system occurs also in between the instants relevant 
to the SM. Given that u/ « 1, this latter variation may 
be considered in adiabatic approximation and it is of the 
order of h [l3], [23[ . As follows from the rough estimates 
above and from the accurate consideration below, the 
variation of energy within the SM dynamics for systems 
of type I is logarithmically larger i.e. larger by the fac- 
tor ln(l//i). Therefore, the variation of energy in be- 
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tween the instants relevant to the SM may be neglected 
at the leading-order approximation for systems of type I 
(the parameter of smallness of the asymptotic theory is 

l/\n(\/h)): AE(-) ~ AE ( s 7n ■ For the sake of notational 
compactness, we shall omit the subscript "sm" further 
in this section. 

For the system (14), the separatrix energy is equal to 
1, while the asymptotic (for E — > E s ) dependence uo{E) 



1(E) ~ I(E S ) 



E R -E 



In 



32 



E. — E 



1 



(17) 



As for the asymptotic value q n (E — > E s ), it is easy to 
see that q(E — > E s ,ip), as a function of ip, asymptotically 
approaches a "square wave", oscillating between — n and 
7r, so that, for sufficiently small j, 



u(E) 
E s = 



In(32/|£7, - E\) 
1, \E. - E\ « 1. 



(15) 



Let us consider the range of energies below E s (the 
range above E s may be considered analogously) and as- 
sume that Uf is close to one of the odd multiples of ojI ■ 
The nonlinear resonance dynamics of the slow variables 
in the range of the approximately resonant energies may 
be described as follows (HH (cf. also US US, 0,0): 



dl dH{I,$) dip dH(I,ip) 



df 



dip 
i 



dt dl 



(16) 



H(I,tp)= / dl (noj — Wf) — nhq n cos(ip) 

JHB.) 



n(E - E s ) -u f (I - I(E S )) - nhq n cos(ip) , 



I ee 1(E) = f 4* 

V ' J E ^u(E)' 
ip = nip — wyt, 



ip = 7r + siga(p)u>(E) 



2tv 



E = H (p,q), 
dq 

liu( E) y/2(E-U(q)) 



2irL 



q n = q n (E) = — J dip q(E,ip) cos(nip), 
|nw-a;/| <w, n=2j-l, j = 1,2,3,..., 

where / and ip are the canonical variables action and 
angle respectively @, S EL S S 13 i ^min is the minimal 
energy over all q,p, E = Hq(p, q); q mnl (E) is the minimal 
coordinate of the conservative motion with a given value 
of energy E; I is the number of right turning points in the 
trajectory [q(r)\ of the conservative motion with energy 
E and given initial state (qo,Po)- 

The resonance Hamiltonian H(I, ip) is obtained from 
the original Hamiltonian H transforming to action-angle 
variables I — ip, with a further multiplication by n; ex- 
tracting the term ujfl (that corresponds to the transfor- 
mation ip — > ip = nip — ui/t); and neglecting all the fast- 
oscillating terms (their effect on the dynamics of slow 
variables is small: see the estimate of the corrections in 
Sec. VI below) i.e. keeping only the resonance term in 
the double Fourier expansion of the perturbation. 

Let us derive the asymptotic expression for 1(E), sub- 
stituting the asymptotic expression (15) for uj(E) into the 
definition of 1(E) (16) and carrying out the integration: 



q^E ^ E s ) ~ (-1) 

Q2j = 0, 
, 7 = 1,2,...« 



2j-l 



(18) 



2u(E) ' 



The next issue is the analysis of the phase space of 
the resonant Hamiltonian (16). Substituting H (16) into 
the equations of motion (16), it is easy to see that their 
stationary points have the following values of the slow 
angle 



V- = o, 



(19) 



while the corresponding action is determined by the equa- 
tion 



nuj - u f T nh -JJ- = °> n = 2 3 - !' ( 20 ) 

where the sign "t" corresponds to ip^i (19). 

As usual (cf. @, S 0, IE S S l24l3i|). the term cx h 
in (20) may be neglected in the leading-order approxima- 
tion, and Eq. (20) reduces to the resonance condition 



(2i-lM^)=« /) 
the lowest-order solution of which is 



E s - E, 



00 



32 exp I — 



(2 j ~ 1)tt 



(21) 



(22) 



Eqs. (19) and (22) together with (17) explicitly deter- 
mine the elliptic and hyperbolic points of the Hamilto- 
nian (16). The hyperbolic point is often called "saddle" 
and corresponds to ip+ or ip- in (19) for even or odd j 
respectively. The saddle point generates the resonance 
separatrix. Using the asymptotic relations (17) and (18), 
we obtain that the resonance Hamiltonian (16) takes the 
following asymptotic value in the saddle: 



TP E 1 ^) 

tisaddle — — ih 

7T 

^^32expf-^^l)-2^. (281 

7T \ Ulf ' 
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The second asymptotic equality in (23) takes into account 
the relation (22). 

As explained in Sec. II above, AE(~\(Of) possesses a 
local maximum at u)f for which the resonance separatrix 
is tangent to the lower GSS curve (Fig. 1(a)). For the 
relevant frequency range (Of — > 0, the separatrix split 
(which represents the maximum deviation of the energy 
along the GSS curve from E s ) approaches the following 
value [J] , in the asymptotic limit h — > 



separatrix (Fig. l(a,b)), which occurs at i\> shifted by 7r 
from the saddle. The condition of equality of H at the 
saddle and at the minimum of the resonance separatrix 
is written as 

Hsaddle = H{I(E ( J) n ), i> saddle + 7T). (30) 
Let us seek its asymptotic solution in the form 



2wh, 



(Of « 1. 



(24) 



As it is shown further down, the variation of energy along 
the relevant resonance trajectories is much larger. There- 
fore, in the leading-order approximation, the GSS curve 
may be simply replaced by the separatrix of the unper- 
turbed system i.e. by the horizontal line E — E s or, 
equivalently, / = I(E S ). Then the tangency occurs at 
ijj shifted from the saddle by n, so that the condition of 
tangency is written as 



'■(idle 



H(I = I(E s ),4> = i> saddle + w)=2h. (25) 



Substituting here H sadd i e (23), we finally obtain the 



U) 



following transcendental equation for u> 



iav{I) = m_A, Is si_pjL. (26) 

^max 

Fig. 3(b) demonstrates the excellent agreement between 
Eq. (26) and the results of simulations for the Hamilto- 
nian system in a wide range of h. 

In the asymptotic limit h — ► 0, the lowest-order explicit 
solution of Eq. (26) is 



E„-E^ n =AE' l 



(l + y)(E s -E^) 



(1 + y)32exp 



7T(2j - 1) 
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(31) 



Substituting (31) and (23) into Eq. (30), we obtain for 
y the following transcendental equation: 



(1 + y) ln(l +y)-y= ——— Xf exp(x f ), (32) 



x f 



where u>m2x is given by Eq. (26). 

Eqs. (31) and (32) describe the left wing of the j-th 
peak of AE('((Jf). Fig. 3(a) demonstrates the good 
agreement between our analytic theory and simulations 
for the Hamiltonian system. 

As follows from Eq. (26), Eq. (32) for = w^L 
reduces to the relation ln(l + y) = 1, i.e. 



7r(2j - 1) 



y>o, 



(33) 



(J) 



(gj - 1)tt 



., = 1.2.... hi(i). !27i 



As follows from Eq. (26), the value of E s - E^ j) (22) 
for io f = WmL is 



Its leading-order expression is: 



4irh 
7W' 

^max 



(28) 



As follows from Eqs. (33), (31) and (28), the maximum 
for a given peak is: 



AeO) = E - E^) (toW ) 



4neh 

7W' 

^max 



(34) 



Fig. 3(c) shows the excellent agreement of this expres- 
sion with the results of simulations for the Hamiltonian 
system in a wide range of h. 

The leading-order expression for AE^lx is: 



2j-l V h 



(29) 



If (Of < then, in the chaotic layer, the largest 

deviation of energy from the separatrix value corresponds 

(i) 

to the minimum energy E„J„ on the nonlinear resonance 

™ mm 



AE {3) 

— lJ - J Tn a\ 



Aeh 



ln(8(2j-l)/h), h^O, (35) 



that confirms the rough estimate (12). 

As (Of decreases, y increases exponentially sharply, as 
follows from Eq. (32). In order to understand how AE^ 
decreases upon decreasing (Of, it is convenient to rewrite 
Eq. (31) expressing the exponent by means of Eq. (32): 
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u f (hi(l + y)-y/(l + y)y 



(36) 



It follows from Eqs. (32) and (36) that AE^ de- 
creases power-like rather than exponentially when w/ is 

decreased. In particular, AE\^ oc l/(wmL — w/) for the 
far part of the wing. 

As for the right wing of the peak, i.e. for Uf > Wmax, 
over the chaotic layer the largest deviation of energy 
from the separatrix value corresponds to the minimum of 
the resonance trajectory tangent to the GSS curve (Fig. 
1(c)). The value of ?p in the minimum coincides with 
tpsaddie- In the leading-order approximation, the GSS 
curve may be replaced by the horizontal line / = I(E S ), 

so that the tangency occurs at tp = tpsaddie + Then 

(1) 

the energy at the minimum E)„ „ can be found from the 
equation 



H{I{E s )^ saddle + tt) = H(I{E^ n ), rjj saddle ) (37) 
Let us seek its asymptotic solution in the form 



e.-eVI = aeP=z(e.-eP) 

n(2j - 1) 



~ 232 exp 
< z < 1, z ~ 1. 



(38) 



Substituting (38) into Eq. (37), we obtain for z the 
following transcendental equation: 



z(l + ln(l/z)) 



8(2j 



- 1) 



— x f exp(xf) (39) 
< z < 1, 



where Wmlx is given by Eq. (26). 

Eqs. (38) and (39) describe the right wing of the j-th 
peak of AE^\ujf). Fig. 3(a) shows the good agreement 
between our analytic theory and simulations. 

As follows from Eq. (26), the solution of Eq. (39) for 
u)f — ► Wmlx is z — ► 1, so the right wing starts from the 
value given by Eq. (28) (or, approximately, by Eq. (29)). 

Expressing the exponent in (38) from (39), we obtain 
the following equation 



AE ( r j) (uf) = 



4nh 



^(l + ln(l/z))- 



(40) 



It follows from Eqs. (39) and (40) that AEf' decreases 
power-like rather than exponentially for increasing uif. 

In particular, AE^ oc l/(wt — uJn-lx) in the far part of 
the wing. 



The further analysis of the asymptotic shape of the 
peak is done in Sec. VII below. 

Beyond the peaks, the function AE^~\u)f) is log- 
arithmically small in comparison with the maxima of 
the peaks. The functions AE^\ujf) and AE^\ujf) in 
the ranges beyond the peaks are also logarithmically 
small. Hence, nearly any combination of the functions 
AEp\u f ) and AE^ j+1) (uj f ) which is close to AE ( r j) {uj f ) 
in the vicinity of Wmlx and to AE^ +1 ^ (ujf) in the vicinity 

of Wmn' may be considered as an approximation of the 
function AE^\ujf) with a logarithmic accuracy with re- 
spect to the maxima of the peaks, AE^lx and ABmi 1 ' , 
in the whole range [wmlx, Wma']. One of the easiest com- 
binations is the following: 



AE^\uj f ) = AE\ 1 \uj f ) for Uf < W « x 
AE { -\u f ) = max{A£«(w/), A£^' +1) (w/)} 



for < uj f < uj ij+1 \ 

^max ^ J ^max ' 



J = 1,2,. 



2^ 



(i) 



(41) 



We used this function in Fig. 3(a), and the analogous 
combination will be also used in the other cases. 

In fact, the theory may be generalized in such a way 
that Eq. (41) would well approximate AE^~\loj) in the 
ranges far beyond the peaks with a logarithmic accu- 
racy even with respect to AE^ [uj f) itself rather than 

to AiJmlx only (cf. the next section). However, we do 
not do this in the present case, being interested primarily 
in the leading-order description of the peaks. 

Finally, we demonstrate in Fig. 4 that the lowest-order 
theory describes quite well the layer boundaries even in 
the Poincare section rather than only in energy/action. 



V. ASYMPTOTIC THEORY FOR SYSTEMS OF 
TYPE II. 

We shall consider two characteristic examples of type 
II, corresponding to the classification given in Sec. III. 
As an example of the system where the separatrix of the 
unperturbed system possesses a single saddle, we shall 
consider the ac-driven Duffing oscillator Mill. As 
an example of the system where the separatrix possesses 
more than one saddle while the perturbation takes equal 
values at the saddles, we shall consider the pendulum 
with an oscillating suspension point @, H, i, HI] ■ The 
treatment of these cases is similar in many respects to 
the one presented in Sec. IV above. So, we present it in 
less details, emphasizing the differences. 



A. AC-driven Duffing oscillator. 

Consider the following archetypal Hamiltonian 0, d, H 
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FIG. 4: A few characteristic Poincare sections in the 2n- 
interval of the energy-angle plane for the system (14) with 
h = 10~ 6 andw/ equal to: (a) 0.236 (maximum), (b) 0.21 (left 
wing), (c) 0.25 (right wing). Results of the numerical integra- 
tion of the equations of motion for the original Hamiltonian 
(14) are shown by red dots. The NR separatrix calculated 
in the leading-order approximation (i.e. by the integration of 
the resonant equations of motion (16) in which lo(E), 1(E) 
and qi(E) are approximated by the explicit formulas (15), 
(17) and (18) respectively) is drawn by the black solid line. 
The NR trajectory (calculated in the leading-order approx- 
imation) tangent to the line E = E s is drawn by the blue 
dashed line. The outer boundary (marked by a thicker line) 
is approximated by: the lower part of the NR separatrix in 
the cases (a) and (b), and by the tangent NR trajectory in the 
case (c) The boundary of the island of stability in the cases 
(a) and (b) is approximated by the tangent NR trajectory 
(which coincides in the case (a) with the NR separatrix). 



H = H + hV, 

2 2 4 

v l q l q 4 

tin = — — + — , 



(42) 

V = — qcos(idft), ft<Cl. 



2 2 ' 4 

The asymptotic dependence of u>(E) on E for E below 



the separatrix energy E s = is the following 0, [H| 

2n 



uj(E) 
E, = 



ln(16/(£ s -£))' 
0, 0<J5, -£<1. 



(43) 



Correspondingly, the resonance values of energies (de- 
termined by the condition analogous to (21)) are 



E,-E' r ''= ... ■i-'^f) 



J = 1,2,3, 



(44) 



The asymptotic dependence of 1(E) is 



1(E) ~ I(E S ) 



E s -E 
2tt 



In 



16 



1 



(45) 



The nonlinear resonance dynamics is described by the 
resonance Hamiltonian H which is identical to Eq. (16) 
in form. Obviously, the actual dependencies u(E) and 
1(E) are given by Eq. (43) and (45) respectively. The 
most important difference is in qj(E): instead of a non- 
zero value (see (18)), it approaches as E — ► E s . Namely, 
it is oc u(E) @,[35l: 



J = 1,2,... « 



uj(EY 



(46) 



i.e. qj is much smaller than in systems of type I (cf. (18)). 
Due to this, the resonance is "weaker" . At the same time, 
the separatrix split S is also smaller, namely ~ hu>f (cf. 
[24l ]) rather than ~ h as for the systems of type I. That is 
why the separatrix chaotic layer is still dominated by the 
resonance dynamics while the matching of the separatrix 
map and nonlinear resonance dynamics is still valid in 
the asymptotic limit h — > 24]. 

Similarly to the previous section, we find the value of 
H in the saddle in the leading-order approximation [36| : 



E s — Er 

2V 



0) 



(47) 



where E s — E^ is given in (44). 

As before, the maximum width of the layer corresponds 
to ojf, for which the resonance separatrix is tangent to 
the GSS curve (Fig. 1(a)). It can be shown [11] that the 
angle of tangency asymptotically approaches ipsaddie + 
7r = 7r while the energy still lies in the resonance range, 
where w(E) w w; ^ m ujf/j. Using the expressions for 
H(E,4>) (cf. (16)), 1(E) (45), qj (E) (46), and taking 
into account that in the tangency E < 6 ~ hujf <C h, the 
value of H at the tangency reads, in the leading-order 
approximation, 



Efangency — f 



v/2' 



(48) 
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Allowing for Eqs. (47) and (48), the condition for the 
maximum, H saddle = H tangency , reduces to 



(49) 



Thus, these values E s —Er are logarithmically smaller 
than the corresponding values (28) for systems of type I. 

The values of u>f corresponding to the maxima of the 
peaks in AE^~\ujf) are readily obtained from (49) and 
(44): 



2nj 



ln(4V2/(7r/i)) : 



j = l,2,...«ln(l/h). (50) 



The derivation of the shape of the peaks for the chaotic 
layer of the separatrix map in the leading order, i.e. 
within the nonlinear resonance (NR) approximation, is 
similar to that for type I. So, we present only the results, 
marking them with the subscript "NR" . 

The left wing of the jth peak of AE^^cof) is described 
by the function 



AE$ H ( W/ ) = 16(l + y)«p (-^ 
2vr\/27i 



ln(l + y)-y/(l + y) : 



-7 



(51) 



where y is the positive solution of the transcendental 
equation 



/ \ i / -i \ n h ( 27TJ 

[l+y)ln{l+y)-y = ^=exp I 



V > 0. (52) 



4V2 V "f J ' 
Similarly to the type I case, 1 + y(iOmix) = e, so that 



As follows from Eqs. (49) and (53), the typical varia- 
tion of energy within the nonlinear resonance dynamics 
(that approximates the separatrix map dynamics) is oc h. 
For the Hamiltonian system, the variation of energy in 
between the discrete instants corresponding to the sep- 
aratrix map @, H, H, 0, 0, HI is also oc h. Therefore, 
unlike the case of type I, one needs to take it into account 
even at the leading-order approximation. Let us consider 
the right well of the Duffing potential (the results for the 
left well are identical), and denote by tk the instant at 
which the energy £ at a given fc-th step of the separatrix 
map is taken: it corresponds to the beginning of the k-th 
pulse of velocity 0, [24j i.e. the corresponding q is close 
to a left turning point qu p in the trajectory [g(r)]. Let us 
also take into account that the relevant frequencies are 
small so that the adiabatic approximation may be used. 
Thus, the change of energy from tk up to a given instant 
t during the following pulse of velocity (t — tk ~ 1) may 
be calculated as 



AE = / dTqhcos(ujfT) ~ hcos(ujftk) / drq 
= h cos(u>ft k ) (q(t) - qup) (56) 

For the motion near the separatrix, the velocity pulse 
corresponds approximately to ip = (see the definition 
of ip (16)). Thus, the corresponding slow angle is ij> = 

jlp - U>ft k ~ -UJftk- 

For the left wing of the peak of AE^\ujf) (includ- 
ing the maximum of the peak too), the boundary of 
the chaotic layer of the separatrix map is formed by the 
lower part of the NR separatrix. The minimum energy 
along this separatrix occurs at i\) = n. Taking this into 
account, and also that "0 — — w e conclude that 

cos(iVftk) — — 1. So, AE < 0, i.e. it does lower the mini- 
mum energy of the layer of the Hamiltonian system. The 
maximum lowering occurs at the right turning point q r t p - 



AE, 



max, NR 



e(£ s -£W(w«J)~2W27i. (53) 



Eq. (53) confirms the rough estimate (13). 

The right wing of the peak is described by the function 



max(|Ai£|) ~ h(q rtp - qu p ) = v2h. 



(57) 



We conclude that the left wing of the j-th peak is de- 
scribed by the following formula: 



AE^ NR {u f ) = 16zexp 
2iry/2h 



_2nf\ 
J 



l + ln(l/z)' 

where z < 1 is the solution of the transcendental equation 



2(1 + ln(l/z)) = -j^= exp (j^j , < z < 1. (55) 
Similarly to the type I case, z(u>f — * Wmlx) — * 1. 



Afif ) ( W/ )~A£g fi ( W/ ) + V2/ l) w f < u&l, (58) 

where AE^fw/) is given by Eqs. (51)- (52). In partic 
ular, the maximum of the peak is: 



AE£l x ~ (27re + l)V2h w 25.6/i. 



(59) 



For the right wing of the peak, the minimum energy 
of the layer of the separatrix map occurs at ip coinciding 
with ip saddle (Fig- l( c )) i- e - equal to 0. As a result, 
cos(u)ftk) — 1 and, hence, AE > 0. So, this variation 
cannot lower the minimal energy of the layer for the main 
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part of the wing, i.e. for LOf < u)^ nd where to^J nd is 
defined by the condition AE^ NR = max(| AE\) ee y/2h. 
For u)f > wH^j, the minimal energy in the layer occurs at 
tp = ir, and it is determined exclusively by the variation 
of energy during the velocity pulse (the NR contribution 
is close to zero at such ip). Thus, we conclude that there 



is a bending of the wing at w/ 



to 



AEV)( Uf ) = AE^^f), < < < nd 

AE^(u f ) = V2h, w/>w 



0) 

bend ' 



UJ 



(i) 

bend 



2nj 



ln(8\/2//i) + 1- 2?r' 



(60) 



where AE^ N r {oj f) is given by Eqs. (54) and (55). 

Analogously to the previous case, AE^~\ojf) may be 
approximated in the whole frequency range by Eq. (41) 
with AE\ j) and AE^ ] given by Eqs. (58) and (60) re- 
spectively. Moreover, unlike the previous case, now the 
theory accurately describes also the range far beyond the 
peaks: AE^ is dominated in this range by the veloc- 
ity pulse contribution AE, which is accurately taken into 
account both by Eq. (58) and by Eq. (60). 

Fig. 5 shows a very reasonable agreement between 
theory and simulations, especially for the 1st peak (37j . 



B. Pendulum with an oscillating suspension point 

Consider the archetypal Hamiltonian 0, H, [H, HH 
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H = H a 



H n = 



P 



hV, 

-cos(g), 



V = — cos(g) cos(u ft), 



h<€ 1. 



(61) 



Though the treatment is similar to the previous case, 
there are also characteristic differences. One of them 
is the following: although the resonance Hamiltonian is 
similar to the Hamiltonian (16), instead of the Fourier 
component of the coordinate, q n , there should be the 
Fourier component of cos(g), V n , which can be shown to 
be: 



FIG. 5: An archetypal example of type II: ac driven Duffing 
oscillator (42). Comparison of theory (solid lines) and sim- 
ulations (circles): (a) the deviation AE^~'(ojf) of the lower 
boundary of the chaotic layer from the separatrix, normalized 
by the perturbation amplitude h, as a function of the pertur- 
bation frequency ujf, for h = 10~ 6 ; the theory is by Eqs. (41), 
(50), (51), (52),' (54), (55), (58) and (60); (b) the frequency of 
the 1st maximum in A£'~'(u/) as a function of h; the theory 
is by Eq. (50); (c) the 1st maximum in AE^~ '{ojf)/h as a 
function of h; the theory is by Eq. (59). 



V 2j ~ {-\y +l -uj{E), E s - E « 1, (62) 

7T 



= 0, 

,? = 1,2,...« 



2vr 
^[EY 



V n ee — / d?/>cos(<j)cos(m/>). 
2tt Jo 

The description of the chaotic layer of the separatrix 
map at the lowest order, i.e. within the NR approxima- 
tion, is similar to that for the ac-driven Duffing oscillator. 



So, we present only the results, marking them with the 
subscript "NR" . 

The frequency of the maximum of a given j-th peak is: 



ln(4//i)' 



j = l,2,...«ln(4// l ). 



(63) 



This expression well agrees with simulations for the 
Hamiltonian system (Fig. 6(b)). To logarithmic accu- 
racy, Eq. (63) coincides with the formula following from 
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(b) 




(c) 



o o o o 



20 - 



FIG. 6: An archetypal example of type II: pendulum with an 
oscillating suspension point (61). Comparison of theory (solid 
lines) and simulations (circles): (a) the deviation AE^~\uj /) 
of the lower boundary of the chaotic layer from the separatrix, 
normalized by the perturbation amplitude h, as a function of 
the perturbation frequency ujf, for h — 10 -6 ; the theory is by 
Eqs. (41), (63), (64), (65), (67), (68), (71) and (73); (b) the 
frequency of the 1st maximum in AE^~'(uf) as a function 
of h; the theory is by Eq. (63); (c) the 1st maximum in 
Aiv - ' (ujf)/h as a function of h; the theory is by Eq. (72). 



Eq. (8) of [18| (reproduced in [23j as Eq. (21)) taken in 

the asymptotic limit h—*0 (or, equivalently, uMLc — > 0). 
However, the numerical factor in the argument of the 
logarithm in the asymptotic formula following from the 
result of [H, is half our value: this is because the 
nonlinear resonance is approximated in (l8l . |23| by the 
conventional pendulum model which is not valid near the 
separatrix (cf. our Sec. Ill above). 

The left wing of the jib. peak of A£^(w^) is described 



by the function 



A25$ fl ( W/ ) = 32(1 + y) exp ( 



Sh 



ln(l+y)-y/(l + y)' 



(64) 



where y is the positive solution of the transcendental 
equation 



(l + y)In(l+»)-y= I ex p(^)' ^ >0 ' ( 65 ) 
Similarly to the previous cases, l+j/(wmL) = e. Hence, 



A 

^^max, NB, ~ e 



(E s -E^(^l))=8eh. 



(66) 



Eq. (66) confirms the rough estimate (13). 

The right wing of the peak is described by the function 



^ ) NR (u f ) = 32zexp(-^j 



8h 



OJf > U) 



(67) 



l + ln(l/z)' 

where z < 1 is the solution of the transcendental equation 



z(i+Hi/z)) = ^ P (^i 

4 V uif 



< z < 1. 



(68) 



Similarly to the previous cases, z(u>f — > Wmax) — > 1- 

Consider now the variation of energy during the veloc- 
ity pulse. Though the final result looks quite similar to 
the case with a single saddle, its derivation has some char- 
acteristic differences, and we present it in detail. Unlike 
the case with a single saddle, the pulse may start close 
either to the left turning point or to the right turning 
point, and the sign of the velocity in such pulses is op- 
posite 0, |24|. As concerns the angle tp in the pulse, it 
is close to — 7r/2 or 7r/2 respectively. So, let us calculate 
the change of energy from the beginning of the pulse, tk, 
till a given instant t within the pulse: 



AE = - 



f drqhdV/dq = h [ drq( 
Jt k Jt k 



— sin(g) cos(oj/t)) 



~ hcos(u)ftk) / drg(— sin(q)) = h cos(uj ftk) cos(q) | \ k 

~ hcos(uj f t k )(cos(q(t)) - 1). (69) 

Here, the third equality assumes adiabaticity while the 
last equality takes into account that the turning points 
are close to the maxima of the potential i.e. close to a 
multiple of 27r (where the cosine is equal to 1). 
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The quantity AE (69) has the maximal absolute value 
at q = 7r. So, we shall lurthcr consider 



VI. ESTIMATE OF THE NEXT-ORDER 
CORRECTIONS 



A£ max = -2hcos(u!ft k ) = -2hcos(2jip k - ij> k ) 



= (-iy +1 2/icos(^) 



(70) 



The last equality takes into account that, as mentioned 
above, the relevant ip k is either — tt/2 or tt/2. 

For the left wing, the value of ip at which the chaotic 
layer of the separatrix map possesses a minimal energy 
corresponds to the minimum of the resonance separatrix. 
It is equal to tt or if the Fourier coefficient V<y is posi- 
tive or negative, i.e. for odd or even j, respectively: see 
Eq. (63). Thus AE'max = — 2h for any j and, therefore, 
it does lower the minimal energy of the boundary. We 
conclude that 



AE®(u,f)~AEfr R ( Uf )+2h, u f (71) 

where AE^ r (uj /) is given by Eqs. (64)-(65). In partic- 
ular, the maximum of the peak is: 



(i) 



AE^>l x ~ (4e + l)2h « 23.7h. 



(72) 



The expression (72) confirms the rough estimate (13) 
and well agrees with simulations (Fig. 6(c)). At the same 
time, it differs from the formula which can be obtained 
from Eq. (10) of [3 (using also Eqs. (1), (3), (8), (9) of 
[III) in the asymptotic limit h — > 0: the latter gives for 
AEiHx the asymptotic value 32h. Though the result [l8[ 
(referred also in [23|]) provides for the correct functional 
dependence on h, it is quantitatively incorrect because 
(i) it is based on the pendulum approximation of the 
nonlinear resonance while this approximation is not valid 
in the vicinity of the separatrix (see the discussion of this 
issue in Sec. Ill above), and (ii) it does not take into 
account the variation of energy during the velocity pulse. 

The right wing, analogously to the case of the Duff- 



ing oscillator, possesses a bending at Uf 



,00 



at 



which AE, 



r,NR 



\AEr, 



2h, that corresponds to 



the switching of the relevant ip by 7r. We conclude that: 



AE^ 



AE, 



U) 



,(J) 



AE^(u f ) = 2h, 
2nj 



(j) 

beiid 



ln(16/ft) 



(73) 



where AE^ NR {ujf) is given by Eqs. (66) and (67). 

Similarly to the previous case, both the peaks and the 
frequency ranges far beyond the peaks are well approxi- 
mated by Eq. (41) with AE^ and AEy' given by Eqs. 
(71) and (73) respectively (Fig. 6(a)). 



We have explicitly calculated only the leading term AE 
in the asymptotic expansion of the chaotic layer width. 
The explicit calculation of the next-order term AE^ next ^ 
is possible but it is rather complicated and cumbersome: 
see the closely related case with two separatrices fzij . 
where most of the next-order contributions are calcu- 
lated quantitatively [38(. In the present paper, where 
the perturbation amplitude h in the numerical examples 
is 4 orders of magnitude smaller than that in [24| , there is 
no particular need to calculate the next-order correction 
C quantitatively. Let us estimate it just qualitatively, 
with the main purpose to demonstrate that its ratio to 
the lowest-order term does vanish in the asymptotic limit 

We shall consider separately the contribution AE^ ext ^ 
stemming from the various corrections within the reso- 
nance approximation (16) and the contribution /\ St E < j: next ^ 
stemming from the corrections to the resonance approx- 
imation. 

The former contribution may be estimated similarly 
to the case considered in [24|: it stems from the devia- 
tion of the GSS curve from the separatrix (this deviation 
reaches 5 at certain angles: see Eq. (7)), from the dif- 
ference between the exact resonance condition (20) and 
the approximate one (21), etc. It can be shown that 

the absolute value of the ratio between AE^ ext ^ and the 
leading term is logarithmically small (cf. [241]): 



\AE, 



(next) I 



AE 



\n(l/h)' 



(74) 



Let us turn to the analysis of the contribution 
^-^{next) ^ . g ^ e CO ntribution stemming from the cor- 
rections to the resonance Hamiltonian (16). It is conve- 
nient to consider separately the cases of the left and right 
wings of the peak. 

As described in Sees. IV and V above, the left wing 
corresponds in the leading-order approximation to for- 
mation of the boundary of the layer by the separatrix of 
the resonance Hamiltonian (16). The resonance approx- 
imation (16) neglects time-periodic terms while the fre- 
quencies of oscillation of these terms greatly exceed the 
frequency of eigenoscillation of the resonance Hamilto- 
nian (16) around its relevant elliptic point i.e. the elliptic 
point inside the area limited by the resonance separatrix. 
As is well known 0, H & H IE @], 

fast-oscillating terms 
acting on a system with a separatrix give rise to the onset 
of an exponentially narrow chaotic layer in place of the 
separatrix. In the present context, this means that the 
correction to the maximal action / stemming from fast- 
oscillating corrections to the resonance Hamiltonian, i.e. 



AE\ 



(next) 



is exponentially small, thus being negligible in 



comparison with the correction A.E-, 



(next) 



see (74)). 
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The right wing, described in Sees. IV and V above, 
corresponds in leading-order approximation to the for- 
mation of the boundary of the layer by the resonance 
trajectory tangent to the GSS curve. For the part of the 
right wing exponentially close in frequency to the fre- 
quency of the maximum, the tangent trajectory is close 
to the resonance separatrix, so that the correction stem- 
ming from fast-oscillating terms is exponentially small, 
similarly to the case of the left wing. As the frequency 
further deviates from the frequency of the maximum, the 
tangent trajectory further deviates from the resonance 



separatrix and the correction AE } 



(next) 



differs from the 



exponentially small correction estimated above. It may 
be estimated in the following way. 

It follows from the second-order approximation of the 
averaging method 39] that the fast-oscillating terms 
lead, in the second-order approximation, to the onset of 
additional terms h 2 Tj(I,ifi) and h 2 T^{I, ip) in the dy- 
namic equations for slow variables I and if) respectively, 
where T/(J, if>) and Tx(I, ip) are of the order of the power- 
law-like function of 1/ ln(l//i) in the relevant range of /. 
The corresponding correction to the width of the chaotic 
layer in energy may be expressed as 



^■^(next) _ 



dt h 2 TjLj(I), 



(75) 



where < m ; n and i max are instants corresponding to the 
minimum and maximum deviation of the tangent trajec- 
tory from the separatrix of the unperturbed system (cf. 



Figs. 1(c) and 4(c)). 
estimated as follows: 



The interval t n 



\<i>>\ 



t min may be 



(76) 



where < ip > is the value of ip averaged over the tangent 
trajectory. It follows from (16) that 



\<i>>\~Wf- lo{E s - 5) 



lo(E s - S) 
\n(l/h) 



LJ 



\n"(l/h) 

(77) 

Taking together Eqs. (75)-(77) and allowing for the 
fact that Tj is of the order of a power-law-like function 
of 1/ ln(l//i), we conclude that 



Thus, we conclude that, both for the left and right 
wings of the peak, (i) the correction AE^ nextS> is deter- 
mined by the correction within the resonance approxima- 
tion AEw' ext \ and (ii) in the asymptotic limit h — * 0, the 
overall next-order correction is negligible in comparison 
with the leading term: 



| A £(nexi)| ^ | A £ 



(next) 



AE 



(next) I 



\AE, 



(next) I 



AE 



AE 



AE 



1 



0. 



(79) 



ln(l/7i) 

This estimate well agrees with results in Figs. 3-6. 

VII. DISCUSSION 

In this section, we briefly discuss the following issues: 
1) the scaled asymptotic shape of the peaks, 2) peaks in 
the range of moderate frequencies, 3) steps in the ampli- 
tude dependence of the layer width, 4) an application to 
the global chaos onset. 

1. Let us analyse the scaled asymptotic shape of the 
peaks. Consider first systems of type I. Then the 
peaks are described in the leading-order approxi- 
mation exclusively within the separatrix map dy- 
namics (approximated, in turn, by the NR dynam- 
ics). As follows from Eqs. (32), (34), (36), (39) 
and (40), most of the peak with a given j can be 
written in the universal scaled form: 



AE^(u) f ) = AE^S f *V / ~ w mL) ) , (80) 



^max J 



where the universal function 5(a) is strongly asym- 
metric: 



Si(a) 



for a<0, 
a) for a>0, 



(81) 



1 



S r (a) 



e(ln(l + y)-y/(l + y))' 
(1 + y) ln(l + y) - y = exp(-a), 
1 



e(l + ln(l/z))' 
z(l + ln(l/z)) = exp(-a). 



AE. 



(next) 



h 2 P(ln(l/h)), 



(78) 



where P{x) is some power-law-like function. 



The value AE. 



(next) 



is still asymptotically smaller than 



the absolute value of the correction within the resonance 
approximation, | AE^ ext ^ |, which is of the order of h or 
/i/ln(l//i) for systems of type I or type II respectively. 



It is not difficult to show that 



Si(a = 0) = 1. 
dSija = 0) = 
da 



1 - e" 



SAa — > — oo) oc : — -, 
\a\ 



S r (a +0) = e-\ (82) 
dS r (a^ +0) 
da 

S r (a — > oo) oc — . 

a 
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Thus, the function S(a) is discontinuous at the 
maximum. To the left of the maximum, the func- 
tion relatively slowly approaches the far part of 
the wing (which falls down power-like) while, to 
the right of the maximum, the function first drops 
jump-wise by a factor e and then sharply ap- 
proaches the far part of the wing (which falls down 
power-like). 

As follows from Eqs. (80), (81), (82) and (27), the 
peaks are logarithmically narrow, i.e. the ratio of 
the half- width of the peak, Aw' J '\ to Wmlx is loga- 
rithmically small: 



~JTT 



ln(8(2j-l)/fc)' 



(83) 



We emphasize that the shape (81) is not restricted 
to the example (14): it is valid for any system of 
type I. 

For systems of type II, the contributions from the 
NR and from the variation of energy during the 
pulse of velocity, as concerns the h dependence, are 
formally of the same order but, numerically, the 
latter contribution is typically much smaller than 
the former one. Thus, typically, the function (81) 
well approximates the properly scaled shape of the 
major part of the peak for systems of type II too. 

2. The quantitative theory presented in the paper re- 
lates only to the peaks of small order n i.e. in the 
range of logarithmically small frequencies. At the 
same time, the magnitude of the peaks is still sig- 
nificant up to the frequencies of the order of one. 
This occurs because, for the motion close to the 
separatrix, the order of magnitude of the Fourier 
coefficients remains the same up to logarithmically 
large numbers n. The shape of the peaks remains 
the same but their magnitude decreases typically 
(but, in some cases, it may even increase in some 
range of frequencies) . The quantitative description 
of this decrease as well as the analysis of more so- 
phisticated cases requires a generalization of our 
theory, that will be presented elsewhere. 

3. Apart from the frequency dependence of the layer 
width, our theory is also relevant for the ampli- 
tude dependence: it describes the jumps [2(| in 
the dependence of the width on h and the transi- 
tion between the jumps and the linear dependence. 

The values of h at which the jumps occur, h^ mp , 
are determined by the same condition which de- 
termines WmL in the frequency dependence of the 
width. The formulas relevant to the left wings of 
the peaks in the frequency dependence describe the 



ranges h > h$ mp 



while the formulas relevant to the 

0) 



right wings describe the ranges h < hj 



ump ' 



4. Finally we note that, apart from systems with a 
separatrix, our work may be relevant to nonlin- 
ear resonances in any system. If the system is 
perturbed by a weak time-periodic perturbation, 
then nonlinear resonances arise and their dynam- 
ics is described by the model of the auxiliary time- 
periodically perturbed pendulum @, H, 0, H, @, 0, Hi • 
If the original perturbation has a single harmonic, 
then the effective perturbation of the auxiliary pen- 
dulum is necessarily a high-frequency one, and 
chaotic layers associated with the resonances are 
exponentially narrow @, H, H, [H, H, 0, [1] while our 
results are irrelevant. But, if either the amplitude 
or the angle of the original perturbation is slowly 
modulated, or if there is an additional harmonic 
of a slightly shifted frequency, then the effective 
perturbation of the auxiliary pendulum is a low- 
frequency one [24] and the layers become much 
wider [40| while our theoretical approach becomes 
relevant. It may allow to find optimal parameters 
of the perturbation for the facilitation of the onset 
of global chaos associated with the overlap in en- 
ergy between different-order nonlinear resonances 
0]: the overlap may be expected to occur at a 
much smaller amplitude of perturbation in com- 
parison with that one required for the overlap in 
case of a single-harmonic perturbation. 



VIII. CONCLUSIONS 

We have further developed the new [24| approach to 
the treatment of the separatrix chaos in the range of 
logarithmically small and moderate frequencies, where 
the chaos takes the largest possible area in phase space. 
The approach is based on the matching between the dis- 
crete chaotic dynamics of the SM and the continuous 
regular-like dynamics of the resonance Hamiltonian. Us- 
ing this approach and taking also into account the dy- 
namics in between instants corresponding to the SM, 
we have presented the first ever accurate asymptotic de- 
scription of high sharp peaks of the width of the sep- 
aratrix chaotic layer in energy as function of the fre- 
quency of a weak time-periodic perturbation, including 
in particular the absolute maximum of the function. Our 
work provides the accurate base to explain former nu- 
merical and heuristic results and intuitive assumptions 
[H, Ho|, [H, [H [!lj2iL corrects the errors of a previous 
heuristic theory 

ism discovers new important fea- 
tures, and opens up new horizons for future studies and 
applications. 

The observed peaks arise due to the involvement of the 
nonlinear resonance dynamics into the separatrix chaotic 
motion. In the context of the magnitude of the peaks, all 
systems are classified into two types: the magnitude of 
the peaks is proportional to the perturbation amplitude 
h times either a logarithmically large factor tx ln(l//i) 
(for systems of type I) or a numerical factor (for systems 
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of type II) . Type I includes systems for which the sepa- 
ratrix of their unperturbed Hamiltonian has more than 
one saddle while the perturbation is not identical on adja- 
cent saddles (an example is an ac-driven pendulum). All 
other systems belong to type II. The latter type includes, 
in particular, a pendulum with an oscillating suspension 
point, for which our result differs from the result [l8j 
since the latter (i) is based on the conventional approx- 
imation of the nonlinear resonance (not valid near the 
separatrix), and (ii) does not take into account the vari- 
ation of energy during the velocity pulse (i.e. in between 
the instants relevant to the separatrix map). Our theory 
is verified by computer simulations. 

The shape of the peaks is strongly asymmetric. In the 
asymptotic limit of small amplitudes, the shape of the 
peaks for type I is universal. For type II, the shape is 
quite similar, differing only by a typically small contri- 
bution stemming from the variation of energy during the 
velocity pulse. 

Our theory describes the jumps of the width as a func- 
tion of the perturbation amplitude h as well as the tran- 



sition between the jumps and the linear dependence. 

Finally, our work suggests a new method for the facil- 
itation of global chaos onset due to the enhanced overlap 
of nonlinear resonances. The theoretical approach devel- 
oped by us may be used to derive the optimal choice of 
parameters of the perturbation leading to the facilitation. 



Acknowledgments 

This work was partly supported by the grant within 
the Convention between the Institute of Semiconduc- 
tor Physics and University of Pisa for 2008 and by the 
Royal Society grant 2007/R2-LJP. We acknowledge dis- 
cussions with Vassili Gelfreich, Igor Khovanov and Oleg 
Yevtushenko. We especially appreciate numerous stim- 
ulating discussions with George Zaslavsky and his role 
for the subject of Hamiltonian chaos on the whole. We 
dedicate this paper to the memory of George Zaslavsky 
who recently passed away. 



[1] G.M. Zaslavsky and N.N. Filonenko, Sov. Phys. JETP 
27, 851 (1968). 

[2] B.V. Chirikov, Phys. Rep. 52, 263 (1979). 

[3] A.J. Lichtenberg and M.A. Lieberman, Regular and 
Stochastic Motion (Springer, New York, 1992). 

[4] G.M. Zaslavsky, R.D. Sagdeev, D.A. Usikov and A. A. 
Chernikov, Weak Chaos and Quasi-Regular Patterns 
(Cambridge University Press, 1991). 

[5] G.M. Zaslavsky, Physics of Chaos in Hamiltonian sys- 
tems (Imperial Colledge Press, 2007). 

[6] G.M. Zaslavsky, Hamiltonian Chaos and Fractional Dy- 
namics (Oxford University Press, 2008). 

[7] S.S. Abdullaev, Construction of Mappings for Hamilto- 
nian Systems and Their Applications (Springer, Berlin, 
Heidelberg, 2006). 

[8] V.G. Gelfreich, V.F. Lazutkin, Russian Math. Surveys 
56, 499 (2001). 

[9] G. N. Piftankin, D.V. Treschev, Russian Math. Surveys 
62, 219 (2007). 

[10] L.E. Reichl, The transition to chaos in conservative clas- 
sical systems: quantum manifestations (Springer, New 
York, 1992). 

[11] Chaos and Stability in Planetary Systems, Lecture Notes 

in Physics, v. 683, eds. R. Dvorak, F. Freistetter, and J. 

Kurths (Springer, Berlin, 2005). 
[12] A.I. Neishtadt, Sov. J. Plasma Phys. 12, 568 (1986). 
[13] Y. Elskens and D.F. Escande, Nonhnearity 4, 615 (1991). 
[14] A.I. Neishtadt, V.V. Sidorenko, and D.V. Treschev, 

Chaos 7, 2 (1997). 
[15] S.M. Soskin, O.M. Yevtushenko, R. Mannella, Phys. Rev. 

Lett. 95, 224101 (2005). 
[16] S.M. Soskin, O.M. Yevtushenko and R. Mannella, 

Commun. Nonlinear Sci. Numer. Simulat., in press; 

doi:10.1016/j.cnsns.2008.06.025. 
[17] V. Gelfreich, private communication. 
[18] LI. Shevchenko, Phys. Scr. 57, 185 (1998). 
[19] A.C.J. Luo, K. Gu, R.P.S. Han, Nonlinear Dyn. 19, 37 



(1999). 

[20] S.M. Soskin, R. Mannella, M. Arrayas and A.N. 
Silchenko, Phys. Rev. E 63, 051111 (2001). 

[21] A.C.J. Luo, Appl. Mech. Rev. 57, 161 (2004). 

[22] V.V. Vecheslavov, Tech. Phys. 49, 521 (2004). 

[23] LI. Shevchenko, Phys. Lett. A 372, 808 (2008). 

[24] S.M. Soskin, R. Mannella, O.M. Yevtushenko, Phys. Rev. 
E 77, 036221 (2008). 

[25] S.M. Soskin, R. Mannella, O.M. Yevtushenko, in Chaos, 
Complexity and Transport: Theory and Applications, 
Proceedings of the CCT-07, eds. C. Chandre, X. 
Leoncini, and G. Zaslavsky (World Scientific, Singapore, 
2008), pp. 119-128. 

[26] V. Rom-Kedar, Physica D 43, 229 (1990). 

[27] Formally, sgn(_E — E s ) is not defined for E = E B but, if to 
shift E from E„ for an infinitesemal value, sgn(_B — E s ) 
acquires the value equal to either +1 or —1, dependently 
on the sign of the shift. Given that Ui+i is proportional 
to sgn(£' 3 — Ei+i) while the latter is random-like (as it 
has been shown above), <7i+i is not correlated with Oi if 
-Ei+i = E s ± 0. 

[28] For the SM relating to the ac-driven spatially periodic 
systems, the time during which the SM undergoes a 
regular-like evolution above the separatrix diverges in the 
adiabatic limit luj — > [l(|, and the width of the part of 
the SM layer above the separatrix diverges too. However, 
we do not consider this case here since it is irrelevant to 
the main subject of the present paper i.e. to to the in- 
volvement of the resonance dynamics into the separatrix 
chaotic motion. 

[29] The GSS curve corresponds to the step of the SM which 
follows the state with E = E„, as described above. 

[30] Namely, that of these two intervals is relevant in which 
the derivative dE/dt in the nonlinear resonance equa- 
tions (see Eq. (16) below) is positive or negative, for the 
case of the upper or lower GSS curve respectively. 

[31] I.e. E determined by Eq. (7) for any t[> except from the 



18 



close vicinity of multiples of n. As shown in [2J|, Eq. (7) 
is irrelevant to the boundary of the chaotic layer in the 
range of ip close to multiples of 7r while the boundary in 
this range of ip still lies in the resonance range of energies, 
where ui(E) fa w^K 

[32] As an example of the heuristic estimate, we use the for- 
mula from 0|: AE^^/h = 2-irujf / cosh(7ru;//2). 

[33] The adiabatic estimate for AE^~'(ujf) is equal approx- 
imately to 27T. The estimate following from the results 
of the mathematical work for w/ ~ 1 is of the same 
order, so that it is conditionally represented in the inset 
in Fig. 2 by the same line as for the adiabatic estimate 
(dashed line). 

[34] S.M. Soskin, R. Mannella and P.V.E. McClintock, Phys. 

Rep. 373, 247 (2003). 
[35] M.I. Dykman, S.M. Soskin, M.A. Krivoglaz, Physica A 

133, 53 (1985). 

[36] The only essential difference consists in that q n at the 
saddle is described by Eq. (46) rather than by Eq. (18). 

[37] The disagrement between the theory and simulations for 
the magnitude of the 2nd peak is about three times larger 
than that for the 1st peak, so that the height of the 
2nd peak is about 30% smaller than that by the asymp- 



totic theory. This occurs because, for energies relevant 
to the 2nd peak, the deviation from the separatrix is 
much higher than that for the 1st peak. Due to the latter, 
the Fourier coefficient q2{E) for the relevant E is signifi- 
cantly smaller than that by the asymptotic formula (42). 
Besides, the velocity pulse contribution AE also signifi- 
cantly decreases while the separatrix split increases as uif 
becomes ~ 1. The generalization of the theory for ujf ~ 1 
will be presented elsewhere. 

[38] The paper [24[ quantitatively calculates all contributions 
into the next-order correction, except the contribution 
stemming from the second-order approximation of the 
averaging method [39]], i.e. from the next-order correction 
to the resonance Hamiltonian. 

[39] N.N. Bogolyubov, Yu. A. Mitropolsky, Asymptotic Meth- 
ods in the Theory of Nonlinear Oscillators (Gordon and 
Breach, New York, 1961). 

[40] This should not be confused with the widening occuring 
with the separatrix chaotic layer in the original pendulum 
if an originally single-harmonic perturbation of a high 
frequency is completed by one more harmonic of a slightly 
shifted frequency: see [221 ] and references therein. 



